Systematic Comparison of Jet Energy-Loss Schemes in a Realistic Hydrodynamic 

Medium 
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We perform a systematic comparison of three different jet energy-loss approaches. These in- 
clude the Armesto-Salgado- Wiedemann scheme based on the approach of Baier-Dokshitzer-Mueller- 
Peigne-Schiff and Zakharov (BDMPS-Z/ASW), the Higher Twist approach (HT) and a scheme 
based on the approach of Arnold-Moore- Yaffe (AMY). In this comparison, an identical medium 
evolution will be utilized for all three approaches: not only does this entail the use of the same 
realistic three-dimensional relativistic fluid dynamics (RFD) simulation, but also includes the use 
of identical initial parton-distribution functions and final fragmentation functions. We are, thus, in 
a unique position, not only to isolate fundamental differences between the various approaches, but 
also to make rigorous calculations for different experimental measurements using "state of the art" 
components. All three approaches are reduced to a version which contains only one free tunable 
parameter, this is then related to the well known transport parameter q. We find that the parame- 
ters of all three calculations can be adjusted to provide a good description of inclusive data on Raa 
versus transverse momentum. However, we do observe slight differences in their predictions for the 
centrality and azimuthal angular dependence of Raa vs. pr- We also note that the value of the 
transport coefficient q in the three approaches to describe the data differs significantly. 



I. INTRODUCTION 

The first seven years of operations at the Relativistic 
Heavy-Ion Collider (RHIC) , performing collisions of gold 
nuclei at Vs NN = 130 GeV and y/s NN — 200 GcV, have 
yielded a vast amount of interesting and sometimes sur- 
prising results [l], 0, 0, [If . Many of these have not yet 
been fully evaluated or understood by theory. There ex- 
ists mounting evidence that RHIC has created a hot and 
dense state of deconfined QCD matter with properties 
similar to that of an ideal fluid @ - this state of matter 
has been termed the strongly interacting Quark- Gluon- 
Plasma (sQGP). 

RHIC has generated a wealth of experimental data on 
high momentum hadron emission, including, but not lim- 
ited to, the nuclear modification factor Raa, its modifica- 
tion as a function of the reaction plane (a measure of the 
azimuthal anisotropy of the cross section) and a whole 
array of high-px hadron- hadron correlations. In these 
observables, one compares the ratio of certain yields in 
a heavy-ion collision to those in a p-p collision, cither 
scaled up by the number of expected binary collisions, 
e.g., for the single hadron suppression factor Raa, or 
directly, as in the case of triggered distributions of asso- 
ciated hadrons, e.g., the Iaa [1, 0,0]- Experimental data 
for most of these observables exist as functions of rapid- 
ity and centrality, for a wide range of pr of the produced 
particle or particles. 

The emission of hadrons with large transverse momen- 
tum is observed to be stro ngly suppressed in central col- 
lisions of heavy nuclei 0, [l(|. The origin of this phe- 



nomenon, commonly referred to as jet-quenching, can 
be understood in the following way: during the early 
pre-equilibrium stage of the relativistic heavy-ion col- 
lision, scattering of partons which leads to the forma- 
tion of deconfined quark-gluon matter often engenders 
large momentum transfers which leads to the forma- 
tion of two back-to-back hard partons. These traverse 
the dense medium, losing energy and finally fragment 
into hadrons which arc observed by the experiments. 
Within the framework of pcrturbativc QCD, the pro- 
cess with largest energy loss of a fast parton is gluon 
radiation induced by collisions with the quasi-thermal 

medium 0, El, El M, EE El, Ei M El ■ 

Computations of jet modification have acquired a cer- 
tain sophistication as regards the incorporation of the 
partonic processes involved. However, the role of the 
medium has often been relegated to the furnishing of an 
overall density and its variation with time [2Cl |2ll 122, |23| . 
Notable departures from these simple treatments in- 
clude attempts to incorporate radial expansion, both 
schematically [13, HH as well as within a fireball evolution 
model [11]. The first attempt to incorporate energy loss 
in a three dimensional (3-D) relativistic fluid dynamical 
(RFD) simulation was carried out by Hirano and Nara in 
Ref. J27[. In this effort, while a full 3-D RFD simulation 
was used, the energy loss of hard jets was carried out 
rather schematically. This approach was also extended 
to the case of two particle correlations in Ref. (HJ. In 
a later effort the authors also incorporated a simplified 
version of the Gyulassy-Levai-Vitev (GLV) energy loss 
formalism at leading order in opacity [29(. 
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Besides the simplified version of the GLV formalism 
used, the authors attempted to apply the results to the 
region in px < 6 GeV, which is the region where data 
were available at the time. In spite of the success of 
Ref. [29l | in explaining the suppression of single inclu- 
sive pions, such a formalism cannot address the flavor 
dependence of the elliptic flow in this region of px- It 
has since been established that jet fragmentation in vac- 
uum is not the primary mechanism of hadronization in 
the range of pr < 6 GeV and there is a sizeable compo- 
nent which arises from recombination. Current rigorous 
implementations of jet modification in dense matter re- 
quire that the pr of the detected hadron be above 6 GeV. 
This allows for a treatment where the final hadronization 
may be treated using the standard vacuum fragmentation 
functions and the ability of a given energy loss formalism 
to compare with experimental data is dependent solely 
on the details of the interaction of the parton with the 
medium in that formalism. This allows for a comparison 
between formalisms where all other components of the 
calculation such as the initial parton distribution, the fi- 
nal fragmentation function as well as the space-time pro- 
file of the medium are identical. This article presents the 
first attempt to perform such a comparison between the 
remaining three formalisms: the BDMPS/ASW, the HT 
and the AMY approach. 

Besides just a comparison between formalisms, this 
paper will simultaneously also apply the different for- 
malisms in comparison to data. A realistic comparison 
with data requires a sophisticated model of the medium. 
The availability of a three-dimensional hydrodynamic 
evolution code [30] allows for a much more detailed study 
of jet interactions in a longitudinally and transversely ex- 
panding medium. The variation of the gluon density in 
such a medium is quite different from that in a simple 
Bjorken expansion. This allows for a step-by-step ap- 
proach to the study of jet-medium interactions. Over the 
past year we have already utilized our evolution model 
to provide the time-evolution of the medium produced 
at RHIC for jet energy-loss calculations performed in the 
BDMPS/ASW [H|, HT [13 and AMY [H approaches. 
In each of the three projects, the inclusive as well as the 
azimuthally differential nuclear suppression factor Raa 
of pions was studied as a function of their transverse mo- 
mentum pt- In addition, the influence of collective flow, 
variations in rapidity, and energy-loss in the hadronic 
phase were addressed for the selected approaches. 

In this manuscript, we shall perform a system- 
atic comparison of jet energy-loss calculations in the 
BDMPS/ASW, HT and AMY approaches. Since we use 
the same medium evolution in all three approaches we are 
in a position to isolate differences among the three calcu- 
lations solely due to their energy-loss schemes. This will 
allow us to answer the question whether the observed dif- 
ferences between the different schemes (when compared 
to data) are due to differing treatment of the medium 
evolution and its coupling to the energy-loss calculation 
or whether they are rooted in more fundamental issues 



related to the energy-loss schemes themselves, e.g. due to 
the approximations and assumptions made when deriv- 
ing the respective schemes. In Sec. II, we briefly review 
the 3D hydrodynamical description of the medium. We 
then discuss in Sec. Ill the theoretical setup of differ- 
ent energy loss schemes and their connection to the 3D 
dynamical evolving medium. Numerical results arc pre- 
sented in comparison to the RHIC data where already 
available in Sec. IV. In Sec. V, we discuss issues related 
to further comparisons of our calculations with the data 
on Raa versus the reaction plane and present concluding 
discussions and an outlook to future work in Sec. VI. 



II. HYDRODYNAMIC DESCRIPTION OF THE 
MEDIUM 

Relativistic Fluid Dynamics (RFD, see e.g. [HI, [H, 
l36l |) is ideally suited for the high-density phase of heavy- 
ion reactions at RHIC, but breaks down in the later, di- 
lute, stages of the reaction when the mean free paths of 
the hadrons become large and flavor degrees of freedom 
are important. The biggest advantage of RFD is that it 
directly incorporates an equation of state as input and 
thus is so far the only dynamical model in which a phase 
transition can explicitly be incorporated. Starting point 
for a RFD calculation is the relativistic hydrodynamic 
equation 

= 0, (1) 

where T^ v is the energy momentum tensor which is given 
by 

T^ v = {e+p)U^U v —pgV. (2) 

Here e, p, U and g^ v arc energy density, pressure, four 
velocity and metric tensor, respectively. The relativis- 
tic hydrodynamic equation Eq. |T]) is solved numerically 
using baryon number ub conservation 

d fl (n B (T,fi)U fl ) = 0. (3) 

as a constraint and closing the resulting set of partial 
differential equations by specifying an equation of state 
(EoS): e = e(p). In the ideal fluid approximation (i.e. ne- 
glecting off-equilibrium effects) and once the initial con- 
ditions for the calculation have been fixed, the EoS is the 
only input to the equations of motion and relates directly 
to properties of the matter under consideration. Ideally, 
either the initial conditions or the EoS should be deter- 
mined beforehand by an ab-initio calculation (e.g. for the 
EoS via a lattice-gauge calculation), in which case a fit 
to the data would allow for the determination of the re- 
maining quantity. Our particular RFD implementation 
utilizes a Lagrangian mesh and light-cone coordinates, 
i.e., (T,x,y,n) where r = y/t 2 — z 2 is the proper time 
and rj is the pseudo-rapidity. This is done in order to op- 
timize the model for the ultra- relativistic regime of heavy 
collisions at RHIC. 
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We assume that hydrodynamic expansion starts at 
t = 0.6 fm. Initial energy density and baryon number 
density are parametrized by 

t(x,y,v) = e max W(x,y;b)H(ri), 
n B (x,y,ri) = n B ma.xW(x,y;b)H(r)), (4) 

where b and e max (fiBmax) are the impact parame- 
ter and the maximum value of energy density (baryon 
number density), respectively. W(x,y;b) is given by 
a combination of wounded nuclear model and binary 
collision model (3?| and H(rf) is given by H(rj) = 
cxp [— (\rj\ — i] n ) 2 /2a 2 ■ 9(\rj\ — 770)] ■ RFD has been very 
successful in describing single soft matter properties at 
RHIC, especially collective flow effects and particle spec- 
tra [30, 1381 KM l40l ] . All parameters of our hydrodynamic 
evolution [30j have been fixed by a fit to the soft sec- 
tor (elliptic flow, pseudo-rapidity distributions and low- 
Pt single particle spectra) , therefore providing us with a 
fully determined medium evolution for the hard probes 
to propagate through. 

III. JET ENERGY-LOSS SCHEMES 

The majority of current approaches to the energy loss 
of light partons may be divided into four major schemes 
often referred to by the names of the original authors: 

• Higher Twist (HT) [M S S, H E HI 

• Path integral approach to the opacity ex pan sion 

(bdmps-z/asw) [13, EE Eil HE El, HE 13, HI 
El, [II 

• Finite temperature field theory approach 

(amy) [II EE dim 

• Reaction Operator approach to the opacity expan- 
sion (GLV) [El [EE [EE m 

All schemes utilize a factorized approach where the fi- 
nal cross section to produce a hadron h with transverse 
momentum px (rapidity between y and y + dy) may be 
expressed as a convolution of initial nuclear structure 
functions (G^(x a ),G^ (ai&), initial state nuclear effects 
such as shadowing and Cronin effect are understood to be 
included) to produce partons with momentum fractions 
x a ,xt, a hard partonic cross section to produce a high 
transverse momentum parton c with a transverse momen- 
tum p and a medium modified fragmentation function for 
the final hadron (D^(z)), 

x — (5) 

dt z 

In the vicinity of mid-rapidity, z = pt/P and t = 
(p ~ x a P) 2 (P is the average incoming momentum of 



a nucleon in nucleus A) . The entire effect of energy loss 
is concentrated in the calculation of the modification to 
the fragmentation function. The four models of energy 
loss are in a sense four schemes to estimate this quan- 
tity from perturbative QCD calculations. While the 
terminology (medium modification) used to describe the 
change in the fragmentation function seems to indicate 
that the medium has influenced the actual process of the 
formation of the final hadrons from the partonic cloud, 
this is not the case. All computations simply describe 
the change in the gluon radiation spectrum from a hard 
parton due to the presence of the medium. The final 
hadronization of the hard parton is always assumed to 
occur in the vacuum after the parton, with degraded en- 
ergy, has escaped from the medium. Note that some 
of the hard gluons radiated from the hard parton will 
also encounter similar "modification" in the medium and 
may endure vacuum hadronization after escaping from 
the medium. Differences between formalisms also arise 
in the inclusion of hadrons from the fragmentation of 
such sub- leading gluons: whereas in approaches which 
compute the change in the distribution of final partons 
(such as AMY) or the change in the distribution of final 
hadrons (such as HT), hadrons from sub-leading gluons 
are implicitly included, formalisms which compute the 
energy loss of the leading parton (such as ASW), do not 
include such sub-leading corrections. 

To better appreciate the approximation schemes, one 
may introduce a set of scales (see Fig. [T|): E or p + , the 
forward energy of the jet, Q 2 , the virtuality of the initial 
jet-parton, p, 1 the momentum scale of the medium and L, 
its spatial extent. Most of the differences between the 
various schemes may be reduced to the different relations 
between these various scales assumed by each scheme as 
well as by how each scheme treats or approximates the 
structure of the medium. In all schemes, the forward en- 
ergy of the jet far exceeds the medium scale, E >> fj,. 
The schemes are presented from one extreme of the ap- 
proximation set (higher twist approach) to the opposite 
extreme (finite temperature approach) , similarities in in- 
termediate steps of the calculation will not be repeated. 
In the following, we shall focus on the first three listed ap- 
proaches, for which we shall present results in section ITVl 
(note that a calculation of GLV jet energy-loss in a 3D hy- 
drodynamic medium has been presented elsewhere [2^|). 



A. Higher Twist Formalism 

The origin of the higher twist (HT) approximation 
scheme lies in the calculations of medium enhanced 
higher twist corrections to the total cross section in 
Deep-Inelastic Scattering (DIS) off large nuclei [58] . One 
re-sums power corrections to the leading twist cross sec- 
tions, which, though suppressed by powers of the hard 
scale Q 2 , are enhanced by the length of the medium. This 
technology of identifying and isolating power corrections 
is used to compute the n-hadron inclusive cross-section. 
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FIG. 1: A schematic picture of the various scales involved in 
the modification of jets in dense matter. 



One assumes the hierarchy of scales E » Q » fi 
and applies this to the computation of multiple Feyn- 
man diagrams such as the one in Fig. [2l this diagram 
represents the process of a hard virtual quark produced 
in a hard collision, which radiates a gluon and then scat- 
ters off a soft medium gluon with transverse momentum 
q_l ~ fx prior to exiting the medium and fragmenting into 
hadrons. At a given order, there exist various other con- 
tributions which involve scattering of the initial quark off 
the soft gluon field prior to radiation as well as scatter- 
ing of the radiated gluon itself. All such contributions 
are combined coherently to calculate the modification to 
the fragmentation function directly. 

The hierarchy of scales allows one to use the collinear 
approximation to factorize the fragmentation function 
and its modification from the hard scattering cross sec- 
tion. Thus, even though such a modified fragmentation 
function is derived in DIS, it may be generalized to the 
kinematics of a heavy-ion collision. Diagrams where the 
outgoing parton scatters off the medium gluons, such 
as those in Fig. [21 produce a medium dependent addi- 
tive contribution to the vacuum fragmentation function, 
which may be expressed as, 



a ^ , 9, dk 2 , a s f 1 dx tt-^ 



In the above equation, Af^j represents the medium 
modified splitting function of parton i into j where a 
momentum fraction x is left in parton j. The argument 
Xl = k 2 _/(2P~p + x{\ — x)) is a momentum fraction de- 
fined such that xlP~ is the formation time of the radi- 
ated parton 1 , where the radiated gluon or quark carries 
away a transverse momentum k±, P~ is the incoming 



momentum of a nucleon in the nucleus and p is the mo- 
mentum of the virtual photon. The scale \x / refers to the 
hard scale of the process. The medium modified splitting 
functions may be expressed as a product of the vacuum 
splitting function Pi_>j and a medium dependent factor, 



APi^ = P^(x) 



C A 2Tra s T^ g (x B ,x L ) 
{kl + {ql))N c mx B y 



(7) 



Where, Ca,N c represent the adjoint Casimir and the 
number of colours. The mean transverse momentum of 
the soft gluons is represented by the factor (qj_). The 
term T^ g represents the quark gluon correlation in the 
nuclear medium, and depends on the four point correla- 
tor, 



(P\m 7 -F-(y 2 )F-' T (y 1 )^(y)\P) 
C{ Pl \^Q) 1 -^{y)\p 1 ){ P2 \F-{y 2 )F-' y {y 1 )\p 2 



(8) 



Where, F~(y 2 ) and F~ a (yi) represent gluon field opera- 
tors at the locations yx, y 2 and ip(y) represents the quark 
field operator. The above correlation function cannot 
be calculated from first principles without making as- 
sumptions regarding the structure of the medium. The 
only assumption made is that that the colour correlation 
length is small. As a result, one may factorize the four 
point function into two separate structure functions, one 
for the original parton produced in the hard scattering 
[this is a quark in Eq. ©] and one for the soft gluon off 
which the parton scatters in the final state. 

While in media with short distance color correlation 
lengths such as the atomic nucleus or a QGP with a large 
Debye mass, this factorization may be generally thought 
to be true, it may fail at very large jet energies where 
saturation effects become important. It should also be 
pointed out that the factorization assumption above falls 
in the same class as the assumption of independent scat- 
tering centers as assumed in the ASW or GLV scheme. In 
the application of this formalism to RHIC data we have 
assumed that the jet energies are not high enough for 
the onset of saturation effects. Another scenario where 
the above factorization may not hold is if there were long 
distance color correlations in the QGP, which have been 
assumed to be absent. If such long distance correlations 
were present then one would have to resort to the defini- 
tion of more general multi-particle operators [such as the 
first line in Eq. |5jl] and parametrize these in comparison 
with experimental data. 

The entire phenomenology of the medium is incorpo- 
rated as a model for the expectation of the second set of 
operators in Eq. (jHJ). This may be characterized in terms 



be referred to using the light cone convention where x^ = (x° ± 
x 3 )/\/2- For the higher-twist scheme, often, x + = (x° + x 3 )/2 
1 Throughout the HT portion of this work, four- vectors will often and x~ = x° — x s . 
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of the well known medium transport coefficient (7(C), at 
location £, where, 



m = 



1 



2tt (2tt) 2 



(9) 



N 2 

1 c 

x exp 

L L V ' 

x (F-i( + e/2,£±/2)F°-(( - £ + /2,-£l/2)) 



o 2 



The Casimir (7r depends on the representation of the 
probe. The transport coefficient is normalized by fitting 
to one data point and a model such as a Woods-Saxon 
distribution for cold matter or 3-D hydrodynamical evo- 
lution for hot nuclear matter is invoked for its variation 
with space-time location. The expectation ( ) is meant to 
be taken in the medium under consideration. Any space 
time dependence is essentially included in the implied 
expectation. 

Closer inspection of Eq. ([9]) reveals that it is a function 
of the jet energy p + . Note that p + is not integrated out. 
The actual dependence on p + depends on the medium 
in question. In the case of confined nuclear media, or 
a quark gluon plasma, the dependence is logarithmic. 
There is also a logarithmic dependence on the virtual- 
ly of the jet which sets in due to radiative corrections 
to the definition in Eq. ([9]). Also, as demonstrated in 
Ref. [1^] , q may even possess a tensorial structure if the 
medium is not isotropic. In the calculations of the cur- 
rent manuscript, both the dependence on the energy and 
virtuality of the jet will be ignored. The medium will be 
assumed to be isotropic. The values of q quoted should 
thus be considered as approximations to the full func- 
tional form. 

Unlike the remaining formalisms, the HT approach is 
set up to directly calculate the medium modified frag- 
mentation function and as a result the final distribution 
of hadrons. This modification to the distribution in- 
cludes both contributions coming solely from the medium 
and those which involve interference between medium in- 
duced and vacuum radiation. The determined constant, 
q, may be used to calculate the average energy loss en- 
countered by a jet. Other advantages of this approach 
include a functional difference between the quark and 
gluon energy loss kernels, i.e., the difference between the 
modification as encountered by a quark jet and a gluon 
jet is not merely assumed to be a ratio of Casimirs, but 
depends strongly on the different splitting and fragmen- 
tation functions. This formalism, offers by far the most 
straightforward generalization to multi-particle correla- 
tions (60| and their modification in the medium. 

A disadvantage of this approach at the current state 
of approximation (similar to the GLV and the ASW but 
different from the AMY approach) is the negligence of 
the quark structure function in the medium: as a result, 
collisions with the medium may not change the flavor 
of the jet parton, however this may continue to occur 
through the splitting kernels. Yet another disadvantage 
is the restriction to single scattering followed by single ra- 




FIG. 2: A typical higher twist contribution used to compute 
the modification of the fragmentation function in medium. 



diation in the medium, which makes this formalism more 
appropriate to thin media. This is partially improved 
by converting Eq. (J7]) to an evolution equation as in (4lj | 
which describes the virtuality evolution of the probe in 
the medium. 

As this formalism is originally cast in cold nuclear mat- 
ter, the applicability of the formalism only depends on 
there being a short distance color correlation length in 
the medium. As a result, it may be used to describe both 
confined and deconfined matter with the inclusion of an 
Ansatz for the variation of q with an intensive property 
of the medium such as energy density e, entropy density 
s or the temperature T and baryon chemical potential 



B. Opacity expansion: Quenching weights 
formalism 



The path integral approach for the energy loss of a 
hard jet propagating in a colored medium was first in- 
troduced in Ref. [131 . ^ was later demonstrated to be 
equivalent to the well known BDMP S ap proach [hi . [l5| 
in the limit of multiple scatterings [6l|. The current, 
most widespread, variant of this approach developed by 
numerous authors [2l], [62j is often referred to as the 
Armesto-Salgado- Wiedemann (ASW) approach. In this 
scheme, one incorporates the effect of multiple scattering 
of the incoming and outgoing partons in terms of a path 
integral over a path ordered Wilson line (47| . 

This formalism assumes a model for the medium as 
an assembly of Debye screened heavy scattering centers 
which are well separated in the sense that the mean free 
path of a jet A ^> 1/fj, the colour screening length of 
the medium [l2j ■ The opacity of the medium n quanti- 
fies the number of scattering centers seen by a jet as it 
passes through the medium, i.e., h = L/X, where L is 
the thickness of the medium. A hard, almost on shell, 
parton traversing such a medium will engender multiple 
transverse scatterings of order fi <C p + . It will, in the pro- 
cess, split into an outgoing parton and a radiated gluon 
which will also scatter multiply in the medium. The 
propagation of the incoming (outgoing) partons as well 
as that of the radiated gluon in this background colour 
field may be expressed in terms of effective Green's func- 
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tions [G(r±, z; r± ', z') (for quark or gluon)] which obey 
the obvious Dyson-Schwingcr equation, 

G(f±,z;r'±,z') = G (f f ±,z;r'±,z') 
- i fdcfd 2 xGo(r±,z;x,QAo(x,OG(x,C;? ± ,z'),(lQ) 



where, Go is the free Green's function and Aq represents 
the color potential of the medium. The solution for the 
above interacting Green's function involves a path or- 
dered Wilson line which follows the potential from the 
location [fj_(z'),z'] to [rj_(z),z]. Expanding the expres- 
sion for the radiation cross section to order A^ n corre- 
sponds to an expansion up to n th order in opacity. 

Taking the high energy limit and the soft radiation 
approximation (x << 1), one focuses on isolating the 
leading behavior in x that arises from the large number 
of interference diagrams at a given order of opacity. As 
a result of the approximations made, one recovers the 
BDMPS condition that the leading behavior in x is con- 
tained solely in gluon re-scattering diagrams. This results 
in the expression for the inclusive energy distribution for 
gluon radiation off an in-medium produced parton as [H| , 

oo oo X X P + 



Co m 



u=r(y) 



dydu 

V=0=r(yi) 



(11) 



where, as always, k± is the transverse momentum of the 
radiated gluon with energy ui and x is a factor that intro- 
duces the kinematic bound. The vectors y and u repre- 
sent the transverse locations of the emission of the gluon 
in the amplitude and the complex conjugate whereas yi 
and yi represent the longitudinal positions. The density 
of scatterers in the medium at location £ is and the 
scattering cross section is a(r). In this form, the opacity 
is obtained as J n(£)d£ over the extent of the medium. 
The Casimir Cr depends on the representation of the jet 
parton. 

Exact analytical expressions for the gluon radiation in- 
tensity distribution are rather involved and only yield 
simple expressions in certain special circumstances. Nu- 
merical implementations of this scheme have focused on 
two separate regimes. In one case, cr(r) is replaced with a 
dipole form qr 2 /n{(,) and one solves the harmonic oscil- 
lator like path integral. This corresponds to the case of 
multiple soft scatterings of the hard probe. In the limit 
of a static medium with a very large length, one obtains 
the simple form for the radiation distribution, 




for it) < lu c , 
for u) > tu r . 



(12) 



Where ui c = J dCC<?(C) is called the characteristic fre- 
quency of the radiation. Up to constant factors, this 
is equal to mean energy lost in the medium ((E)) i.e., 
lu c ~ 2(E) /(a. s Cr). In the other extreme, one expands 
the exponent as a series in na; keeping only the lead- 
ing order term corresponds to the picture of gluon radi- 
ation associated with a single scattering. In this case, 
the gluon emission intensity distribution has been found 
to be rather similar, once scaled with the characteristic 
frequency appropriate for this situation. For dynamical 
medium of finite extent, the characteristic frequency and 
the overall mean transverse momentum gained by the jet 
(qL) will have to be estimated based on an Ansatz for 
the space time distribution of the transport parameter 

q- 

Due to the soft limit i.e., to — > used, multiple gluon 
emissions are required for a substantial amount of energy 
loss. Each such emission at a given opacity is assumed in- 
dependent and a probabilistic scheme is set up, wherein, 
the jet loses an energy fraction AE in n tries with a 
Poisson distribution 12111 . 



P n (AE) 



-(i) 



-IT 



du), 



dl 

' dwi 



S(AE-J2u)ill3) 



where, (I) is the mean number of gluons radiated per co- 
herent interaction set. Summing over n gives the proba- 
bility P(AE) for an incident jet to lose a momentum frac- 
tion AE due to its passage through the medium. This 
is then used to model a medium modified fragmentation 
function, by shifting the energy fraction available to pro- 
duce a hadron (as well as accounting for the phase space 
available after energy loss), 



D(z,Q 2 



i D 
dAEP(AE) — 



(l-AB'^ 2 ) 

1 - AE 



(14) 



The above, modified fragmentation function is then used 
in a factorized formalism as in Eq. ((5|) to calculate the 
final hadronic spectrum. 

In marked contrast to other approaches, this scheme 
presents the advantage of easy interpolation between the 
cases of few hard scatterings and multiple soft scatterings 
and is thus applicable to both thin and thick media. The 
inclusion of the zero opacity term makes this the only 
formalism, to date, which includes interference between 
vacuum radiation and radiation induced by multiple soft 
scattering in the medium. It suffers from the disad- 
vantage of having approximated the medium in terms 
of heavy static scattering centers. As a result, clastic 
energy loss is vanishing in this scheme. As the formal- 
ism is setup to calculate the energy loss probability of 
the leading hard parton, estimation of the change in the 
distribution of final associated (sub-leading) hadrons or 
partons is not straightforward. 

Along with the HT formalism, this approach also ne- 
glects any flavor changing scatterings in the medium. 
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as. 




FIG. 3: A ladder diagram evaluated in the AMY formalism 



Also similar with HT is the treatment of both confined 
and deconfined matter on the same footing: one essen- 
tially makes an Ansatz for the variation of q with an 
intensive variable of the medium e.g., e, s,T, /is- 



C. Finite temperature field theory formalism 

In this scheme, often referred to as the Arnold-Moore- 
Yaffc (AMY) approach, the energy loss of hard jets is con- 
sidered in an extended medium in equilibrium at asymp- 
totically high temperature T — > oo. Owing to asymptotic 
freedom the coupling constant g — > at such high tem- 
peratures, and a power counting scheme emerges from 
the ability to identify a hierarchy of parametrically sepa- 
rated scales T >> gT >> g 2 T etc. In this limit, it then 
becomes possible to construct an effective field theory of 
soft modes, i.e., p ~ gT by summing contributions from 
hard loops with p ~ T, into effective propagators and 
vertices |64j |. 

One assumes a hard on-shell parton, with energy sev- 
eral times that of the temperature, traversing such a 
medium, undergoing soft scatterings with momentum 
transfers ~ gT off other hard partons in the medium. 
Such soft scatterings induce collinear radiation from the 
parton, with a transverse momentum of the order of 
gT. The formation time for such collinear radiation 
~ 1/ (g 2 T) is of the same order of magnitude as the mean 
free time between soft scatterings [5lf . As a result, multi- 
ple scatterings of the incoming (outgoing) parton and the 
radiated gluon need to be considered to get the leading 
order gluon radiation rate. One essentially calculates the 
imaginary parts of infinite order ladder diagrams such as 
those shown in Fig. [3J this is done by means of integral 
equations [65j . 

The imaginary parts of such ladder diagrams yield the 
1 — > 2 decay rates of a hard parton into a radiated gluon 
and another parton. These decay rates are then used 
to evolve hard quark and gluon distributions from the 
initial hard collisions, when they are formed, to the time 
when they exit the medium, by means of a set of coupled 
Fokkcr-Planck like equations for quarks, anti-quarks and 
gluons [22|, HH, H3 , which may be written schematically 



dPj{p,t) 
dt 



ab 

-Pj( P ,t) 



Pa(p+k,t) 

dTj b { P ,k,t) 
dkdt 



dr- b (p + k, P ,t) 



dkdt 



(15) 



In the above equation, j = q,q,g, and we sum over all 
relevant partonic processes for each evolution equation. 
In contrast to all other schemes, this approach also in- 
cludes the absorption of thermal gluons as well as quark 
anti-quark pair annihilation and creation. 

The initial jet distributions are taken from a factorized 
hard scattering cross section as in Eq. ((3J) . In the limit of 
single scattering, these rates may be taken directly from 
the corresponding Gunion-Bertsch cross sections [66| for 
an on-shell parton to radiate a gluon on soft scattering 
with another in-medium parton. 

The final hadron spectrum at high pt is obtained by 
the fragmentation of jets in the vacuum after their pass- 
ing through the medium. In this approach, one calcu- 
lates the medium modified fragmentation function by 
convoluting the vacuum fragmentation functions with the 
hard parton distributions, at exit, to produce the final 
hadronic spectrum 53], 



DHz,r ± ,< 



E 



dp^D^Pip^p^A).^) 



where the sum over j' is the sum over all parton species. 
The two momentum fractions are z = Ph/Pj and z' = 
Ph/Pj> , where pj and pji are the momenta of the hard par- 
tons immediately after the hard scattering and prior to 
exit from the medium, and ph is the final hadron momen- 
tum. P(j>ji\pj, rj_, 0) represents the solution to Eq. (|15|) . 
which is the probability of obtaining a given parton j' 
with momentum py when the initial condition is a parton 
j with momentum pj. The above integral depends im- 
plicitly on the path taken by the parton and the medium 
profile along that path, which in turn depends on the 
location of the origin r± of the jet, its propagation an- 
gle <f> with respect to the reaction plane. Therefore, one 
must convolve the above expression over all transverse 
positions r± and directions (f>. 

The use of an effective theory for the description of 
the medium and the propagation of the jet, makes this 
approach considerably more systematic than the two pre- 
vious approaches: both the properties of the jet and the 
medium are described using the same hierarchy of scales. 
It remains the only approach to date which naturally 
includes partonic feedback from the medium, i.e., pro- 
cesses where a thermal quark or gluon may be absorbed 
by the hard jet 2 . In contrast to ASW and HT, this ap- 



2 While an attempt to include such effects in the higher twist for- 
malism have been made in Ref . |67| , these remain as phenomeno- 
logical extcntions and have not been included in this manuscript. 
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proach also includes flavor changing interactions in the 
medium. Elastic energy loss may also be incorporated 
within the same basic formalism [68|. Note that AMY 
assumes a thermalized partonic medium and neglects the 
quenching of jets in the confined sector. In addition, in- 
terference between medium and vacuum radiations is not 
yet considered. 

The use of HTL effective theory to describe both the 
jet propagation in the medium and the properties of the 
medium itself does suffer from one caveat: this scheme 
approximates the bulk structure of the medium as a 
weakly coupled plasma of quarks and gluons. The per- 
turbativc estimates of the energy density (e) differs from 
the e(T) obtained from lattice calculations (at 3T C > T > 
T c ). The r)/s required to reproduce the observed magni- 
tude of elliptic flow in viscous fluid dynamical simulations 
is at least a factor of two lower than perturbativc results. 
The application of such a scheme to the modification of 
hard jets involves an aspect of phenomenology where the 
coupling constant is used fit parameter. 



D. Geometry and Discussion of the different 
Schemes 



T via T 3 or the entropy density s of the local medium. 
In the analysis presented in this paper we maintain this 
methodology, however, some surprising results of scaling 
the ASW q with T 3 and the HT q with e 3 / 4 will also be 
presented. 

The scaling of q with e 3 / 4 , T 3 or s will by construc- 
tion yield identical results for a QGP with an ideal gas 
equation of state: e = 3p. However, for a more realistic 
non-ideal equation of state as used in our hydrodynamic 
calculation, the value of q will be affected by the choice 
of scaling variable, in particular if energy-loss persists to 
temperatures below Tc- Figure 2J investigates the devi- 
ations from the ideal gas scaling by plotting the normal- 
ized time evolution of temperature T 3 , energy-density 
e 3 / 4 and entropy density s. As can be seen after 2 fm/c 
the curves start to deviate from the ideal gas power law 
behavior and start to show differences for times later than 
t = 3 fm/c. The first order phase-transition contained 
in our equation of state results in a striking difference 
between the temperature and the energy- or entropy- 
density scaling at the critical temperature and below. 
We note that the proper scaling law for q is a priori not 
known, even though we see no reason why it should not 
be calculable in QCD. 



As mentioned previously, all parameters of our hydro- 
dynamic evolution (30l | have been fixed by a fit to the 
soft sector (elliptic flow, pseudo-rapidity distributions 
and 1ow-j»t single particle spectra), therefore providing 
us with a fully determined medium evolution for the 
hard probes to propagate through. The hydrodynamic 
calculation provides a time-evolution of the tempera- 
ture, energy-density, flow velocity and QGP to hadron 
gas fraction within all hydrodynamic cells composing the 
medium through which the hard probes evolve. The in- 
corporation of this information within the different jet 
energy-loss schemes is described in the following subsec- 
tions. 

The mean impact parameters for the different evolu- 
tion sets have been set to b =2.4, 4.5, 6.3, 7.5 fm, cor- 
responding to 0-5, 10-15, 15-20, 20-30 % centrality, re- 
spectively. These values were estimated via the number 
of nucleon-nucleon binary collisions and the number of 
participant nucleons in Ref. In this work, the fo- 
cus will lie on the two extreme centrality bins in the list 
above: the 0-5% bin and the 20-30% bin. All of the RFD 
calculations utilized here have an initial thermalization 
time of To = 0.6fm/c. Any values of parameters, such 
as <7, which are dependent on the bulk properties of the 
medium will be quoted at this time. 

All three energy loss schemes are sensitive to certain 
bulk properties of the evolving matter: while in the case 
of the AMY formalism this is decidedly the temperature, 
the relation between q and the intensive variables of the 
medium in the HT and ASW formalisms is unspecified. 
Traditionally, the q in the ASW scheme has been related 
to the the energy density e via e 3 / 4 while the q in the 
HT scheme has been scaled either with the temperature 
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FIG. 4: (Color online) Time evolution of temperature T 3 , 
energy-density e 3 ^ 4 and entropy density s as a function of 
time t in the central cell of the hydrodynamic evolution for 
Au+Au collisions at RHIC. All curves are normalized to their 
maximum values at r = 0.6 fm/c. 



1. Higher Twist 

In the preceding section, the medium modification of 
the final fragmentation function in the Higher Twist for- 
malism was shown to be dependent on the transport co- 
efficient q [see Eq. ([5])]. In the evolving system formed in 
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the collision of two nuclei, this transport coefficient has 
both a space and time dependence [i.e., q (x, y, z, r)]. 
Phcnomenologically, this dependence is taken to scale 
with some intensive variable of the medium, in this case, 
the dimensionally equivalent quantities of T 3 or the en- 
tropy density s, i.e., 



»/ x - r y±{x,y,z,T)T 3 (x,y,z,T) 
q{x,y,z,T) = q -3 (17) 

1 o 



x [R(x, y, z, t) + c HG {1 - R(x, y, z, r)} ] 



where, T{x,y,z,T), j±(x,y, z,t) and R{x,y,z,r) repre- 
sent the temperature, flow transverse to the jet and the 
volume fraction in the plasma phase at the space-time 
point x,y,z,T. ft is this information that is extracted 
from the RFD simulation. Though the RFD simulations 
start at r = 0.6 fm/c, the values of T and j± at r = 0.6 
fm/c are extrapolated as constants to t = which repre- 
sents the time of the initial hard scatterings (the effect of 
different extrapolation schemes involving linearly rising 
or dropping values of q as r — > has been found to be 
rather small and will not be studied in this effort). The 
factors qo,To represent the maximum q and temperature 
achieved in the simulation; in this particular version of 
RFD, To = 0.405 GeV and go is a fit parameter adjusted 
to fit one data point of the Raa, at one centrality. 

The factor chg may be interpreted in two ways. In 
essence, it accounts for the fact that the quenching in 
the hadronic phase may not be as effective as that in the 
partonic phase at the same temperature. Since the en- 
tropy density in a given phase is proportional to T 3 with 
the constant of proportionality demonstrating a weak de- 
pendence on temperature, chg m ay be tuned to convert 
the scaling of q with T 3 into a scaling with s. This is ap- 
proximately achieved with a chg ~ 0.2, as can be seen in 
Fig. [U Unless specified otherwise, this is the value used 
for chg m all the plots in this paper. Thus chg is not 
a fit parameter and is not tuned to fit any experimental 
data point. It has only three possible choices of chg = 
which corresponds to no quenching in the hadronic phase, 
chg = 1 which corresponds to exact scaling of q with T 3 
and chg = 0.2 which corresponds to approximate scaling 
of q with s. 

Given a choice of chg and the overall fit parameter 
go, we use Eq. ([7]) to calculate a medium modified frag- 
mentation function; then Eq. ([5]) is used to compute the 
total cross section and the nuclear modification factor 
Raa- The overall fit parameter go is tuned to fit one ex- 
perimental data point, at one centrality and px- For the 
current effort, the fit parameter is set by requiring that 
the Raa at px = 10 GeV in the most central event (0-5% 
centrality) is 0.2. With the value of go and chg fixed, the 
variation of Raa as a function of px (integrated over or 
with respect to the angle with the reaction plane) and 
centrality of the collision are predictions. 



2. ASW 

As in the previous case, we have to formulate the en- 
ergy loss problem for a dynamical medium in which the 
transport coefficient q acquires a space and time depen- 
dence. As done in previous calculations within the ASW 
formalism, we use a scaling with the local energy density 
e 3 / 4 along the path £ = (x(r), y(r), z(r), t) of a parton 
as 



q(0 = K-2-e^ 4 (0 



(18) 



This scaling of q is assumed to be valid in both the 
partonic and the hadronic phase. The precise form of the 
path £ is determined once the hard initial vertex (xq, yo) 
in the transverse plane, the outgoing parton rapidity 77 
and the angle of the parton with the reaction plane (f> is 
specified. The parameter K in Eq. (Tj"8|) is regarded as a 
parameter to account for the uncertainty in the selection 
of a s and possible non-perturbative effects increasing the 
quenching power of the medium (see discussion in [69l|). 

Given this spacctime dependence of the transport co- 
efficient along a parton trajectory, the energy loss prob- 
ability distribution can be computed from the two line 
integrals 



/•OO /"OO 

: (r o ,0)=/ and (gX)(r , = / 

Jo Jo 

(19) 



Here, uj c is the characteristic gluon frequency, setting 
the scale of the energy loss probability distribution (see 
expression [T2"l) and (qL) is a measure of the path- length, 
weighted by the local quenching power. Analogously to 
the overall fit parameter go in the HT case, the parameter 
K is fit one data point of the Raa, at one centrality. 

For times prior to r = 0.6 fm/c, i.e. the starting point 
of the RFD simulation, we neglect any medium effects, 
i.e. assume q = 0. Note that for a purely radiative energy 
loss model where the average energy loss grows quadrat- 
ically with pathlength in a constant medium the effect 
of initial time dynamics is systematically suppressed and 
no strong dependence of the energy loss on variations of 
the initial time is observed. 

Using a dynamical scaling law [49| , u> c and (g) can then 
be mapped onto a static equivalent scenario. Using the 
relation R = 2uj^/(qL) as an input, we then determine 
P(AE) using the numerical results from 2l| and com- 
pute the medium-modified fragmentation function from 
Eq. (fT4)) . The resulting expression (valid for a single 
path) must then be averaged over the whole geometry 
with a weight corresponding to the probability of finding 
an initial hard vertex at xo,yo. 
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3. AMY 

The strength of the transition rate in pQCD is con- 
trolled by the strong coupling constant a s (T), tempera- 
ture T and the flow parameter (3 (the velocity of thermal 
medium) relative to the jet's path. The value for the 
coupling constant used (along with the assumption of 
a thermalized partonic medium) may be related to the 
transport coefficient q as derived from computations in 
HT or ASW by either a direct computation of the oper- 
ator product in Eq. or a computation of the mean 
transverse momentum squared per unit length as gained 
by a jet which propagates through the medium without 
radiation. 

In a 3D expanding medium, the transition rate is first 
evaluated in the local frame of the thermal medium, then 
boosted into the laboratory frame, 



dT(p,k,t) 



dkdt 



1 



■0) 



dT(p ,k ,t ) 



lab 



dkodto 



(20) 



local 



where k = fc(l - Vj ■ (3 2 and to = ty/l - /3 2 

arc momentum and the proper time in the local frame. 
As jets propagate in the medium, the temperature and 
the flow parameters depend on the time and the po- 
sitions of jets, and the 3D hydrodynamical calculation 
[30l ] is utilized to determine the temperature and flow 
profiles. The energy loss mechanism is applied at time 
To = 0.6 fm/c, when the medium reaches thermal equi- 
librium, and switched off when the medium reaches the 
hadronic phase. 



IV. APPLICATION TO RHIC DATA 

In the preceding sections, a description of the theoret- 
ical setup underlying each of the three schemes as well as 
the phenomenological connection between them and the 
RFD simulations was expounded upon. In this section, 
we present the results of our numerical calculations. The 
primary quantity of interest will be the nuclear modifi- 
cation factor (Raa) defined as 



Raa = 



d<J AA (b m , n ,b max ) 
dyd 2 p T 



d?bT AA {b) 

d 2 bdyd 2 pT 



d<jPP( PT ,y) 
dyd 2 p T 



(21) 



where, Taa represents the nuclear overlap function which 
is proportional to the number of binary collisions at the 
mean impact parameter (b) . The mean impact parameter 
for a given range of centrality is essentially set by the 
RFD simulation used to calculate the soft observables. 
The Raa is calculated both integrated as defined above 
or as function of the angle with respect to the reaction 
plane. 



The range of pr of the detected hadron is set high 
enough for the applicability of pQCD. In this paper, the 
lower bound is set at pr = 6 GeV. This choice is essen- 
tially dictated by the regime where recombination (70l | 
begins to contribute to the yield. The upper limit is 
set at pt = 20 GeV which represents the highest pr for 
which experimental data exist. The focus in this paper 
will essentially be on two different centralities of (b) = 2.4 
fm which represents the rather central collisions with a 
centrality in the range from — 6% and a somewhat more 
peripheral event with a (b) = 7.5 fm which corresponds 
to the 20 — 30% range of centrality. 
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FIG. 5: (Color online) Nuclear modification factor Raa in 
Au-Au collisions at 0-5% (top) and 20-30% (bottom) central- 
ity calculated in the ASW, HT and AMY approaches com- 
pared to data from PHENIX 

Figure [5] shows the nuclear modification factor Raa as 
a function of px in Au-Au collisions at 0-5% (top) and 20- 
30% (bottom) centrality calculated in the ASW, HT and 
AMY approaches compared to data from PHENIX pTll ] . 
The parameters for the respective calculations are fixed 
to one data point in the 0-5% centrality calculation - 
the dependence on pt and centrality of the nuclear col- 
lision are then predictions by the respective energy-loss 
calculations. As can be seen, the parameters for all three 
approaches (initial maximal value for the transport co- 
efficient go or coupling constant a s in the AMY case) 
can be adjusted such that the approaches are able to 
describe the centrality dependence of the nuclear modi- 
fication factor reasonably well. For a gluon jet, the val- 
ues are qo « 4.3 GeV 2 /fm for the HT approach, qo ~ 
18.5 GcV 2 /fm for the ASW formalism and a s w 0.33 
for the AMY approach, which can be converted into a 
value of % 4.1 GeV 2 /fm. While values of qo, have 
been presented up to the first decimal point, one should 
note that the error involved is never less than the exper- 
imental error (See Sec. V for further discussion on this 
issue). Beyond this, there remain the usual uncertain- 
ties related to using a leading order hard scattering cross 
section, e.g., the choice of the appropriate scale for the 
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structure and fragmentation functions. There also exist 
additional sources of error in the estimations of q which 
arise from the set of approximations used in each of the 
formalisms to reduce the functional dependence on the 
properties of the medium down to one parameter. 

The reader will note a somewhat smaller value of go 
quoted for the HT formalism in Fig. [5] Since, the HT 
approach was originally developed for DIS on a large nu- 
cleus, it has become customary to quote the value of go 
for a quark which is always the produced hard parton 
in DIS (see Refs. HHHfl). Besides this difference, there 
remain various caveats associated with this value of q 
which have been discussed in Sec. Ill [in particular see 
the discussion surrounding Eq. JHJ)] . 

For the case of the ASW formalism, we have used the 



relationship [721 ] . 



So 



2Ke 



3/4 



(22) 



to convert the parameter K in the ASW approach to go- 
In a previous estimate using this formalism [69( , the value 
of qo was quoted to be somewhat lower. This is simply 
due to the earlier time To = 0.6 fm/c at which go is being 
quoted in the current manuscript. In Ref. [69j, To was set 
to 1 fm/c. 

In AMY, the relationship between q and the coupling 
a s reads 



<1 



C A g 2 Tm 2 D ^ g^* 



2tt 



mo 



(23) 



where g™ ax is the largest transverse momentum rele- 
vant for the collinear emission. One estimate is that 
(g™ ax ) 2 « ET, where E is the energy of the jet, and T 
the temperature. Evaluating the above expression for 3 
quark flavors, a s = 0.33, a temperature of 0.4 GeV and 
a jet energy of 20 GeV, one obtains q — 4.1 GeV 2 /fm. 
Even though this formulation is only logarithmic in the 
jet energy, it is however more suggestive than precise 
[Z3. Note that the ASW value for q at r = 0.6 fm/c 
and eo = 55 GeV/fm 3 lies a factor of 3.6 higher than the 
Baier estimate for an ideal QGP, q « 2 • e 3 / 4 [7^|, while 
the AMY estimate is in line with that from Baier, and 
the HT calculation lies about a factor of 1.6 below that 
value. 

The large difference in go values between HT, AMY 
and ASW has been pointed out previously. However, our 
calculation shows for the first time that this difference is 
not due to a different treatment of the medium or ini- 
tial state. Note that the numbers quoted here reflect 
the different medium scaling laws referred to as being 
the natural choices for the respective approaches, namely 
temperature scaling for AMY, energy-density scaling for 
ASW and entropy-density scaling for HT, as discussed 
in the previous section. If we choose to perform the 
jet energy-loss calculation with temperature ~ T 3 scal- 
ing for all three approaches, we find values for go be- 
ing 10 GcV 2 /fm for ASW, 2.3 GcV 2 /fm for HT and 
4.1 GeV 2 /fm for AMY. Likewise, if we employ energy- 
density scaling ~ e 3 / 4 , we find go = 18.5 GeV 2 /fm for 
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ASW 


HT 


AMY 
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go 


go 


go 


T(r» 


10 GeV 2 /mi 


2.3 GeV 2 /fm 


4.1 GeV 2 /fm 




18.5 GeV 2 /fm 


4.5 GeV 2 /fm 




s(r, t) 
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TABLE I: Values of go, i-e., the q at r = To = 0.6 fm/c in 
the cell at r = of the — 5% centrality event, in the dif- 
ferent energy loss schemes. Also presented is the variation 
of go with different choices of scaling of q(r, r) with different 
local intensive properties of the medium; where T(f, t) is the 
temperature, e(f, t) is the energy density and s(f, r) is the en- 
tropy density at location (r, r) in the RFD simulation. Given 
the model of the medium in AMY, q may only be calculated 
as a function of T (see Eq. I23[). hence the entries correspond- 
ing to e and s scaling are left blank. Calculations in the ASW 
scheme with q scaled with s have not yet been performed and 
so the entry for s scaling has been left blank. 



ASW and g = 4.5 GeV 2 /fm for HT (the AMY calcu- 
lation can only be performed utilizing temperature scal- 
ing). Both ASW and HT consistently show a rise of a 
factor of two in go when switching from temperature scal- 
ing to energy-density scaling. The different values for go 
in the different schemes with different choices of scaling 
with T, s and e 3 / 4 are presented in Table. |TJ 
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FIG. 6: (Color online) Raa as a function of azimuthal angle at 
pr = 10 GeV/c and pr = 15 GeV/c for all three approaches 
in the 20-30% centrality bin. 

We find that slight differences appear between the ap- 
proaches when Raa is studied as a function of azimuthal 
angle. This can be seen in Fig. O where Raa is plot- 
ted as a function of azimuthal angle at pt = 10 GeV/c 
and pr = 15 GeV/c for all three approaches in the 20- 
30% centrality bin. Figure[?l shows the same calculation, 
but with all curves normalized by their respective az- 
imuthally averaged Raa ~ we observe that for the pr bins 
chosen, the AMY and HT calculations exhibit the same 
peak-to-valley ratio and shape, whereas the ASW calcu- 
lation shows a more pronounced difference between in- 
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plane and out-of-plane emission. The azimuthal spread 
is insensitive to variation of the transverse momentum, 
which is manifest in the comparison between the solid 
(pr = 10 GeV/c) and the dashed (pr = 15 GeV/c) lines. 
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FIG. 7: (Color online) Raa as a function of azimuthal angle at 
pr = 10 GeV/c and pr = 15 GeV/c for all three approaches 
in the 20-30% centrality bin, normalized by the azimuthally 
averaged value of Raa for the respective calculations 



In order to further quantify the difference between the 
three approaches we calculate the ratio of the out of plane 
Raa over the in plane Raa as a function of transverse 
momentum - this is shown in Fig. [5J We find that AMY 
and HT exhibit the same peak to valley ratio through- 
out the entire range of transverse momenta, even though 
the absolute values for Raa differ by approximately 10%. 
The ASW calculation systematically shows a stronger az- 
imuthal dependence than the HT and AMY calculations 
- the cause of which will require a more detailed analysis 
to determine. 
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FIG. 8: (Color online) Ratio Raa for out of plane vs. in plane 
emission as a function of pr at b=7.5 fm impact parameter 
for all three approaches. 



Note, however, that the agreement in the peak-to- 
valley ratio for AMY and HT does not translate into 
these approaches being identical in terms of the in-plane 
and out-of-plane Raa values vs. pr- Fig. [5] shows that 
systematic differences on the order of 15% exist between 
all three approaches in the absolute value of Raa at fixed 
azimuthal angle as a function of px- 
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FIG. 9: (Color online) Raa for out of plane vs. in plane 
emission as a function of pr at b=7.5 fm impact parameter 
for all three approaches. 

In order to investigate the spatial response of the jet 
energy-loss schemes to the medium, we define the follow- 
ing quantity, 



P(x,y) 



Tab{x,v) ■ Raa{x,v) 
J dxdyT AB (x, y) ■ Raa(x, y) ' 



(24) 



where the local position-dependent nuclear suppression 
factor RAA{x,y) is weighted with the nuclear over- 
lap probability function Tab{x, V)- Figure [T0l shows 
P(x, y = 0) as a function of x for a quenched jet mov- 
ing in the positive x direction through the center of the 
medium (y = 0). 

Integrating the quantity P(x, y) over y yields the es- 
cape probability of a hadron with a transverse momen- 
tum between 6 and 8 GeV/c originating from a quenched 
jet moving in the positive x direction in the transverse 
plane as a function of of its production vertex along the 
x-axis: 



P(x) = / dyP(x,y) 



(25) 



The result is shown in Fig. [II]- it is remarkable how well 
the three different approaches agree with each other in 
this quantity. Since the same hard scattering probabil- 
ity was used as input in all three cases, the agreement 
in P(x) really shows that all three approaches yield the 
same suppression factor as a function of production ver- 
tex of the hard probe, i.e. that they probe the density of 
the medium in the same way. 
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FIG. 10: (Color online) Survival probability P(x,y) of a 
hadron with 6-8 GeV/c or 12-14 GeV/c transverse momentum 
moving along the positive rr-axis at through the center of the 
medium (y = 0) in the transverse plane as a function of x. 
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FIG. 11: (Color online) Escape probability of a hadron with 
6-8 GeV/c transverse momentum moving along the positive 
a;-axis in the transverse plane as a function of x. 



limit of the applicability of jet quenching calculations, 
a comparison may not be completely out of place, but 
would carry large uncertainties with it. 

However, the data from run-2 of the PHENIX col- 
laboration [74| which was used to deduce the Raa 
versus the reaction plane demonstrates an integrated 
Raa of 0.41 ± 0.03(stat) ± 0.06(sys) in the 20-30% 
ccntrality events, in noticeable contrast to the value of 
0.35 ± - OM(stat) ± - OM(sys) as seen in Fig. 
from the run-4 data set. While the two data sets agree 
within systematic errors, the discrepancy between the 
two is too large for a meaningful comparison of our cal- 
culation, which was fit to the run-4 data set. 

An estimate of the variation of the fit parameters re- 
quired to encompass both data sets leads to differences 
of the order of 20 - 40% in q. Plotted in Fig. [H are 
the predictions for the Raa versus reaction plane for the 
standard values of the fit parameters obtained from the 
comparison with the run-4 data set in Fig. [Sj Also plot- 
ted are readjusted plots for the Raa versus the reaction 
plane where the single fit parameters in each of the mod- 
els was tuned such that the integrated Raa in the 20- 
30% centrality bin achieved a value of 0.41. The new 
values of the fit parameters (included in the figure) are 
q = 1.6 GcV 2 /fm for the HT, K = 2.4 for the ASW and 
a s = 0.27 for the AMY calculations respectively. One 
should note that a simple renormalization of our Raa vs. 
4> curves to the data would not be appropriate, since the 
value of q affects the magnitude of the azimuthal spread. 

A detailed and meaningful theory-experiment compar- 
ison, encompassing different data sets as well as their 
respective statistical and systematic errors in a proper 
fashion will require a sophisticated statistical analysis be- 
yond the scope and aim of the work presented here. Such 
an analysis has been demonstrated for one particular the- 
ory calculation of inclusive Raa vs. px compared to one 
experimental data set in [76] . The feasibility of extending 
such an analysis to multiple data sets, observables and 
theory calculations has yet to be determined. 



V. NORMALIZATION AND FURTHER 
COMPARISON TO DATA 

As we have seen in the previous section, there do ex- 
ist noticeable differences in the Raa as a function of the 
azimuthal angle between the three approaches. A com- 
parison to experimental data for this particular observ- 
able would thus constitute an important experimental 
input and possibly serve as a discriminator. Recently, 
data for Raa versus the reaction plane have become 
available in the px = 5 — 8 GeV region [74|. Unfor- 
tunately this px range, which in terms of the data will 
be dominated by the lower px boundary, still sits in the 
region in which particle production is significantly influ- 
enced by parton recombination as hadronization mecha- 
nism |7fJ, l75j . Since we regard px = 6 GeV as the lower 



VI. SUMMARY AND CONCLUSION 

In summary, we have calculated the modification of 
hard jets in a 3D hydrodynamic medium, in three dif- 
ferent approaches which were constrained to use the 
same initial structure functions, the same final vac- 
uum fragmentation functions, the same nuclear geometry 
and identical three dimensional evolution of the produce 
dense matter. In this first, unified, attempt to under- 
stand jet modification in dense matter, the focus was 
restricted to single inclusive observables. The nuclear 
modification factor [Eq. (|2~lj) ] was computed as a func- 
tion of the transverse momentum, centrality of collision, 
as well as, the angle with respect to the reaction plane. 
This was followed by a more detailed, though purely the- 
oretical, analysis of jet origin distribution for the Raa 
as a function of the reaction plane, as well as, the Raa 
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FIG. 12: (Color online) Raa vs. reaction plane in the 20- 
30% centrality event at pr = 6 GeV for different choices of 
the single fit parameters q, K, a 3 

for jet origins restricted to lie on a narrow belt on the 
reaction plane. 

In the comparisons above, both the HT and the ASW 
schemes have been simplified to the point that all predic- 
tions depend on only one tunable parameter: this is the 
(FF) correlator in the HT approach and the K param- 
eter in the ASW approach. In the most rigorous formu- 
lation of AMY, there exist no free parameters except for 
the temperature; this however, has already been spec- 
ified by the RFD simulation. In the phenomenological 
application of the AMY approach used here, the strong 
coupling constant is treated as a parameter: It has, thus, 
been disassociated from the temperature. 

These single free parameters from all three approaches 
were tuned to fit one data point, usually chosen as the 
integrated Raa at 8 GeV in the 0-5% centrality events. 
The data used for this comparison as shown in Fig. O were 
taken from the PHENIX run-4 data set [7l|. Our com- 
parison shows that under identical conditions (i.e. same 
medium evolution, same choice of parton distribution 
functions, scale etc.) all three jet energy-loss schemes 
yield very similar results. This finding is very encour- 
aging since it indicates that the technical aspects of the 
formalisms are well under control. However, we need to 
point out that there still exists a puzzle regarding the ex- 
tracted value for the transport coefficient <?o , which spans 
a factor of four from a value of 2.3 GeV 2 /fm for the HT 
approach on the lower end, to 4.1 GeV 2 /fm for AMY and 
10 GeV 2 /fm for ASW on the high end, when using the 
same temperature scaling law for all three approaches. 
While the discrepancy among these approaches is not 
new, our work has been able to decisively rule out dif- 



ferences in the medium evolution or initial setup as a 
cause for the differing values of q. We are led to conclude 
that these remaining differences are due, to the different 
approximations applied, to the different energy scales in- 
volved, to the different assumptions on the structure of 
the QCD matter, inherent in these different approaches. 

There exist multiple future directions for the system- 
atic and unified approach to jet modification in dense 
matter presented here. Due to the assumption of a ther- 
malized plasma, elastic energy loss may be straightfor- 
wardly included in AMY. Including elastic energy loss, 
however, represents a significant extension to the HT and 
ASW approaches which has, only recently, been under- 
taken and thus, this topic has not been included in the 
comparisons presented in this article. The current ef- 
fort was restricted to single inclusive observables; hence, 
the simplest extension will be to apply a similar analysis 
to both single and multi-particle observables in tandem. 
Such comparisons will undoubtedly lead to stronger con- 
straints on the formalism and hence deeper insights in the 
nature of the theory of jet modification used. Another 
direction is to use a somewhat different initial condition 
and equation of state for the medium evolution. A natu- 
ral extension in this direction is to the study of jet mod- 
ification in viscous fluid dynamical simulations. Viscous 
simulations, necessarily seem to require an initial state 
with greater spatial anisotropy. We believe, that it is in 
this direction that measurements and theoretical calcula- 
tions of the Raa versus the reaction plane will have most 
relevance, as a means to discriminate between different 
initial state profiles. The approximations which have re- 
sulted in the reduction of formalisms such as the HT and 
the ASW to a dependence on only one parameter will 
eventually have to be relaxed. The different parameters 
in these schemes represent actual physical properties of 
the produced matter which may indeed be measurable 
given a detailed and extensive set of experimental mea- 
surements. 
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